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Abstract. - The relaxation time of weakly interacting classical spins is calculated by intro- 
ducing the averages of the local dipolar field, obtained by thermodynamic perturbation theory, 
in a rigorous expression for the single-spin thermoactivation rate in a weak but arbitrarily ori- 
ented field. At low temperatures the non-trivial dependence of the superparamagnetic blocking 
on the damping coefficient, numerically found by Berkov and Corn, is reproduced by our model 
and interpreted in terms of the deviations from uniaxial anisotropy associated to the transversal 
component of the dipolar field acting on each spin. 



Introduction. - The study of single-domain magnetic particles has been an active field 
of research since the pioneer work of Stoner and Wohlfarth who studied the hysteretic 
rotation of the magnetisation over the magnetic-anisotropy energy barrier under the influence 
of applied fields, and Ncel j^] who predicted that at non-zero temperature the magnetisation 
can surmount the energy barrier as a result of thermal agitation. Important progress has 
been made since Brown ||^ derived the Fokker-Plank equation for the probability distribution 
of spin orientations, starting from the stochastic Landau-Lifshitz equation, and calculated 
the relaxation time for uniaxial particles in a longitudinal field. Recent work on spins with 
non-axially symmetric potentials revealed Q] a large dependence of the relaxation time on 
the damping coefficient A in the medium-to-weak damping regime (A measures the relative 
importance of the precession and the damping in the dynamics). Experiments on individual 
nanoparticles analysed with accurate asymptotes of the relaxation time , gave damping 
coefficients in that regime: A « 0.05-0.5. 

Turning our attention from independent particles to systems of interacting particles, the 
complexity of the problem increases drastically, as it becomes a many-body problem with 
a long-ranged and reduced symmetry interaction mechanism — the dipolar interaction. The 
approach to study the relaxation time r has so far been based on how the energy barriers of 
the spins are modified by the interactions 0,1). With the r so obtained, one can study the 
effects of the interaction on the superparamagnetic blocking (the maximum in the dynamical 
response at the temperature where t becomes of the order of the observation time). However, 
the barrier-based approach to determine t corresponds to assuming A — > c», so that dynamical 
features such as the precession of the spins are disregarded, as incisively noted by Berkov and 
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Gorn 1^. Indeed, numerical integration of the stochastic Landau-Lifshitz equation for weakly 
interacting systems revealed Q non-trivial effects of finite damping on the superparamag- 
netic blocking, such as enhanced shifts of the temperature of the maximum of the dynamical 
response and non-monotonic behaviour of its height with the interaction strength. 

In this article we apply thermodynamic perturbation theory (treating the anisotropy energy 
exactly and the dipolar energy perturbatively) to calculate the averages of the dipolar field 
produced at the position of a given spin. These averages are then introduced in a rigorous 
low-field expansion of the single-spin relaxation rate, which can be combined with a recently 
derived perturbative formula for the equilibrium susceptibility [ p^ to provide a model for the 
dynamical susceptibility of weakly interacting spins. In the low temperature range, where 
the superparamagnetic blocking takes place, our model recreates the non-trivial damping 
dependence of the blocking found in Ref. ||^. The analytical treatment permits us to ascribe 
the observed features to the deviations from uniaxial symmetry of the potential due to the 
transversal components of the dipolar field, deviations which render the relaxation time very 
sensitive to the damping. 

Low- field relaxation rate. - In weak fields, we can handle the field dependence of the 
relaxation time t by expanding the relaxation rate F = 1/r in powers of the field components 
(for convenience we use the field in temperature units ^ — inB / k^T , where m is the magnetic 
moment). As the spins have inversion symmetry in the absence of the field, the relaxation 
rate will not change if the field is reversed (since T accounts for jumps over the energy barrier 
in both directions). Therefore, T should be an even function of ^, and hence all odd order 
coefficients in the expansion must be identically zero. For spins with uniaxial anisotropy, 
the second order corrections can be written in terms of the longitudinal ^|| and transverse 
components of the field with respect to the anisotropy axis. Using, for instance, the invariance 
of the relaxation rate upon field reflection through the barrier plane, where does not change 
but is reversed, we see that no quadratic term of the form ^||^^ can occur. This leads to 
the following generic form for the relaxation rate of uniaxial spins in the presence of a weak, 
arbitrarily directed, field 

r^Fo(i-i-c||e| + ciei) , (1) 

where Fq is the zero-field relaxation rate and the expansion is actually valid to third order. 
The coefficients C|| and cj_ will be determined by choosing special configurations in which they 
are known (strictly longitudinal and transverse fields). 

We shall restrict our attention to low temperatures where the superparamagnetic blocking 
takes place. Then, the zero-field relaxation rate is given by Brown's result 

Fo = -^a^/'e-" , (2) 

where td (oc 1/A) is the relaxation time of isotropic spins, a — A/k^T is the anisotropy 
barrier in temperature units, and corrections of order 1/a are disregarded due to the low T 
assumption. Note that in this case (and also in a longitudinal field) the damping parameter 
A only enters through td and hence it only matters to establish a global time scale. In other 
words, the results for different A presented in units of td show complete dynamical scaling, and 
in this sense the A dependence is said to be trivial. Expanding the corresponding expression 
for F in the presence of a longitudinal field ^|| (disregarding again 1/a terms), we find 



F(e||,a =0)~Fo(l + i^f 



p. E. JoNssoN AND J. L. Garci' a-Palacios : Relaxation time of weakly interacting superparamagnets 3 



Comparison with the general expansion gives the longitudinal coefficient C|| — 1/2. 

There is no general expression for the relaxation time in the presence of a non-zero trans- 
verse field valid for all values of the relevant parameters. Nevertheless, Garanin et al. Q have 
derived a low temperature formula valid for weak transverse fields, which is perfectly suited 
for our purpose of determining c^, namely 



1 1 
2^' 2^ 



Here a — Xa^l"^ and 7(0, z) — dtt°'^^ e^* is the incomplete gamma function. Comparing 
with the expansion (|l|) one gets the transverse coefficient c± = F/4. 

On gathering these results we finally get the desired expression for the low temperature 
relaxation rate in weak fields 



i + 5efi + H«Ki , (3) 



which constitutes a straightforward generalisation of the formula of Garanin et al. [Q to an 
arbitrary field orientation. The function F takes into account without further approximations 
the effects of the precession. F decreases towards 1 for strong damping (where C|| and cj^ are of 
the same order of magnitude), while F grows as 1/A for weak damping, where the relaxation 
time turns to be very sensitive to the damping (or to transverse fields). 

Let us compare the rigorous Eq. (^) with the corresponding result of M0rup and co-workers 
(see, for instance, Eq. (6) in Ref. Q), which for I/ct <C 1 can be written as 

rM0rup-ro(i-^ii + ie^ + iei) • (4) 

This expression has the undesirable feature of a term linear in the field, which can be attributed 
to accounting for the escape from one of the potential wells only. Actually, if we average over 
both wells the linear term cancels out and the corrected formula equals the overdamped A — *■ 00 
Hmit of Eq. (||) (then F ^ I). This is natural, since in Ref. Q the effects of the field were 
only considered via barrier changes, so that no gyromagnetic effects were included, and their 
result can only be correct when the precession can be neglected (A — > 00). Besides, as their 
linear term disappeared upon random anisotropy averaging, their results will correspond to 
the limits overdamped plus random anisotropy of those derived here from Eq. (^). 

Averages 0/ the dipolar field. - Let us consider a system of magnetoanisotropic spins 
coupled via the dipole-dipole interaction. The dipolar field at the position ri of the spin Si 
created by all other spins is given by 

where the term with j — i is omitted, m is the magnitude of the magnetic moment, and in 
the dipolar tensor r-ij ~ fi — rj and f ij — rij/rij. All lengths are measured in units of the 
characteristic length a, which is defined in such a way that is the mean volume occupied 
around each spin. In a simple cubic arrangement a is the lattice constant and for nanoparticles 
of volume V the volume concentration of particles is V/a^. For notational simplicity we are 
assuming that the parameters characterising the different spins are identical (it is immediate 
to generalise the expressions for different anisotropy constants, magnetic moments, volumes, 
etc.). 
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For spins with uniaxial anisotropy the total energy of the system can be written as 

- (iE = a ^(si • fiiY + ^ Si • Gi-j ■ Sj , (6) 



i>j 



where (3 — l/k^T, fii is the unit vector along the anisotropy axis of the ith spin, a is the 
anisotropy barrier divided by the thermal energy, and 

The quantity hd is a convenient temperature independent measure of the interaction strength 
(oc concentration) equal to the magnitude of the field, measured in units of the maximum 
anisotropy field Bk — 2A/m, produced at a given position by a spin located at a distance a. 

Using thermodynamic perturbation theory (Ref. |]l2|, § 32) to expand the Boltzmann 
distribution W = exp{—l3E) in powers of ^d, one gets an expression of the form 

Wis,) = W^is,) [1 + edii'i (s.) + •••]. (8) 

Here = Yii '^a.^ exp(— /3i?a) is the Boltzmann distribution of the non interacting en- 
semble (in our case Eg_ includes the magnetic anisotropy) and w\{si) is linear in the dipo- 
lar energy (and hence quadratic in the spins). The calculation of the observables reduces 
to computing averages weighted by the non interacting probability distribution = 
J ^ exp(— /3i?a)(- • •) of low grade powers of the spin variables. In this way, the results 
are exact in the magnetic anisotropy and only perturbational in the dipolar interaction. 

The averages of terms linear and quadratic in s, weighted by the noninteracting distribu- 
tion, can be calculated by means of the following "algorithms" [|l^ 

((s • vi))^ = 0, ((s • vi)is- V2))^ - 1(1 - ^2) VI-V2 + S2 (n ■ vi){n ■ V2) , (9) 

where vi and V2 are arbitrary constant vectors. The first average vanishes (actually any odd 
power) since the magnetic anisotropy has inversion symmetry [Wa(— Si) = Wa(si)] and there is 
no external bias field. The quantity 6*2 is the average of the second order Legendre polynomial 
S2 = {^[3(s ■ — l])a, and it can be written in terms of the one-spin partition function Za, 
as 

S2{o)-li^-^^-l^ Za = V^erf(iV^) . (10) 



2 VcrZa 2(7 J 2 

In the equations for the relaxation rate the field enters squared [Eq. (||)]. Therefore, we 
do not calculate the statistical mechanical average of Bi [Eq. (^)] and plug it into F (which 
would be a sort of mean-field approach) but we average instead the combinations of field 
variables as they enter in the expression for F. Since B\ = — i3y , we average, by means 
of the algorithms (^), the square of the field and the square of the projection along the local 
anisotropy axis B"^.. — {Bi ■ flif , getting (^i = mBi/kBT) 



(^Iw ) = f E [ (1 - ^2) (n. • G,, . G., . ,tO + 3^2 (n. • G., • rt,)^] , (11) 



3 

t2 



3 
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- (1 - ^2) (fii ■ Gij ■ ■ fii) - 3S2 {fii ■ Gjj • fij)"^] . (12) 
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The averages have been calculated to order (so we just need Eq. (||) to zero order; wi only 
enters in the third order corrections). Note that the averaged fields depend on T via S'2(cr), 
reflecting the fact that the field created at a given position can be different if, for example, 
the source spins are almost freely rotating (high T, 52 0) or are almost parallel to their 
anisotropy axes (low T, 52 — > 1). 

One may wonder about the validity of these field averages below the superparamagnetic 
blocking, where the spins are not in complete equilibrium. However, since at those tempera- 
tures the spins are still in quasi-equilibrium confined to one of the two wells, we can repeat 
the derivation of the algorithms (^ restricting the phase space for integration to one well. 
In this case, averages of the form (s • vi)^ do not vanish, and should be considered together 
with ((s - Vi){s- V2)) which being even in s is not modified. The extra terms associated with 
{s ■ vi)g^, however, vanish if the overall state is demagnetised, and we recover Eqs. ([lT|)-(|l2|). 
This can be interpreted as if two nearby blocked spins, one in its upper well and the other in 
the lower, create a net field at the position of a third spin similar to that created by any of 
them when in equilibrium, since each of the spins approximately compensates for the restricted 
phase space of the other. 

The general expressions for the longitudinal and transversal fields are notably simplified in 
some important situations. For a system with parallel anisotropy axes (e.g., in a single crystal 
of magnetic molecular clusters, or ferrofluids frozen in a strong magnetic field) we equate all 
the ftj to n. For a system with randomly distributed anisotropy axes we replace expressions 
involving /(rij) by integrals J (Pn f{n) = /, and use {n ■ vi){n ■ V2) = \vi-V2- In both cases the 
final expressions involve some rapidly convergent sums over the lattice. Let us concentrate in 
cases of "sufficiently isotropic" lattices, in the sense of fulfiUing '^{rx)^ — J2i^y)'^ — X)(''^)'^' 
e.g., cubic and completely random lattices (incidentally, the type of arrangements for which 
in the classical Lorentz cavity field calculation the contribution of the dipoles inside the 
"small" sphere vanishes). Similarly we consider large enough systems, so that all spins have 
approximately equivalent surroundings (then the index i on the different quantities can be 
dropped). Under these circumstances the "lattice sums" involved are 

n^2j2r;^\ T = ^(n.G.,.n)2, (13) 

3 3 

where the terms with j ^ i are of course omitted. 

For a system with aligned anisotropy axes the averaged fields are given in terms of the 
lattice sums (^) by the compact expressions 

{^') = f - ^2) n + 352 T] , (ei) = I [(2 + 52)7^ - 352 T] , (14) 
while for randomly distributed anisotropy axes they read 

'W) = jT^^ (Ci) = |27^. (15) 

In both cases (^2) = + is simply given by the temperature independent result (^2) = 
^^7?.. This independence also holds for the field components after the random anisotropy 
average [Eq. (|l^)], but not for parallel axes [Eq. (p^j where the field components remain T 
dependent. 

Superparamagnetic blocking. - The dependence of the features of the superparamagnetic 
blocking on the interaction strength — ^d/2cr has been a subject of some controversy. Two 
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main approaches addressed this problem on the basis of the modifications of the energy 
barriers by the interactions. However, this type of approach overlooks the fact that not only 
the energy landscape is important, but also how the spin evolves in it, depending on if the 
spin is precessing almost freely (weak damping) or strongly damped. Indeed, Berkov and 
Gorn have shown with rigorous Langevin dynamics simulations that for weak interactions 
the position of the blocking temperature Th (where x" reaches its maximum) of strongly 
damped spins is hardly affected by the interaction strength, whereas for weak damping Tb 
significantly decreases with /id and the peak height behaves non-monotonically. The apparent 
discrepancy with some experimental results, in which increases with h^, was attributed 
in Ref. to the different behaviour of systems with weak anisotropy (or moderate-to-high 
interactions), where the energy barriers are mostly due to the interactions and hence grow 
with hd- 

In order to assess the features of the superparamagnetic blocking emerging from our model, 
we make use of a rigorous perturbative expansion of the equilibrium susceptibility Xoq in 
powers of (to second order, see Ref. pM). This expression, together with the relaxation 
rate T obtained when the averaged fields (|l4|)-(|l5|) are introduced in Eq. (|^), can be combined 
in a simple Debye-type formula 

r 

Naturally, at low enough temperatures | [lO| the results will become invalid by the very na- 
ture of the approximations involved, although this does not affect the characteristics of the 
superparamagnetic blocking for weak interactions. 

The dynamical response for a large spherical sample with parallel anisotropy axes and 
simple cubic lattice structure is shown in Fig. |l| (the lattice sums required are TZ = 16.8 and 
T — 13.4). In the overdamped case, Tb is not noticeably affected by the dipolar interaction 
and the height of the susceptibility peak decreases monotonically with h^- This corresponds to 
the effect of a slight decrease of Tb found in Ref. |^ for random anisotropy, which is observable 
only at very high frequencies (e.g., Mossbauer). For weak damping, however. Fig. |l| shows that 
the blocking temperature significantly decreases as the interaction strength increases and, in 
addition, the peak height of x" initially rises for small values of hd and then decreases for 
larger values. The same behaviour can be seen in x' with the appropriate choice of parameters. 

The features shown in Fig. ^are in complete agreement with the results obtained by Berkov 
and Gorn The analytical treatment employed here, however, allows us to readily trace back 
the origin of the results obtained and to interpret them in terms of the different hd dependences 
of r and Xon- For overdamped systems the coefficients C|| =1/2 and c± — F/4 ~ 1/4 in the 
expression (g) for the relaxation rate are of order unity and lead to a slight increase of F 
(decrease of Tb) with hd, while the entire curve is lowered by the reduction of Xoq with hd- 
For weak damping, however, F oc 1/A is large and makes F very sensitive to the interaction 
strength, moving the peak quickly towards low T as soon as hd departs from zero. Then, as the 
coefficients in Xcq are not as large as c±, the initial decrease of the equilibrium susceptibility 
is smaller than the increase associated to the quick shift of the blocking to low T, where Xcq is 
higher (roughly cx 1/T), producing the rise of the peak. If hd is further increased the decrease 
of Xeq starts to be competitive and the peaks are reduced as they shift to lower temperatures. 

The different behaviours are related to the presence of transverse components of the local 
fields, which create a saddle point in the uniaxial potential barrier of the spins, turning the 
relaxation rate sensitive to A We can picture the underlying physical mechanism as 
follows [|l^. Consider one spin that after a "favourable" sequence of fluctuations, reaches a 
point close to the top of the barrier but does not surmount it. In the subsequent spiralling 
down back to the bottom of the potential well, a strongly damped spin descends almost 
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immediately, whereas a weakly damped spin executes several rotations (~ 1/A) about the 
anisotropy axis. This allows the latter spin to pass close to the saddle area, where it will have 
additional opportunities, not available for the damped spin, to cross the barrier. As we see, 
this mechanism duly combines the characteristics of the potential and the dynamical evolution 
of the spins in the potential. Note finally that the transverse components are non-zero even 
for parallel anisotropy axes [Eq. 0], so there is no need, as in the case of non-interacting 
particles, to appeal to obhque fields, applied Q or probing [Q, to find a large sensitivity to 
the damping in interacting systems. 

Summary and Conclusions. - We have proposed a model for the relaxation time of 
weakly interacting superparamagnets. The single-spin relaxation time at low fields is gen- 
eralised for an arbitrary directed field, and the components of the local field calculated by 
thermodynamical perturbation theory. The non-monotonic behaviour of the height of the 
dynamical susceptibility peak with the interaction strength and the increased lowering of the 
blocking temperature with decreasing damping, discovered numerically by Berkov and Gorn, 
are captured by our model. These features are interpreted in terms of the different sensitiv- 
ity, depending on the damping strength, of the relaxation rate to transverse fields, which for 
interacting spins are provided by the dipolar interaction. 

* * * 
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Fig. 1 - Imaginary component of the dynamical susceptibility vs. temperature (the real component is 
shown in the inset) for a spherical sample and spins placed in a simple cubic lattice. The anisotropy 
axes are all parallel and the response probed along their common direction. The dipolar interaction 
strength hd = (.d/^a is: hd = Q (thick lines), 0.004, 0.008, 0.012, and 0.016 from (a) up to down and 
(b) right to left. The frequency is lotd/ct = 27r x 0.003. 



